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ABSTRACT 


A  complete  analysis  and  explicit  solution  is  presented  for  the  problem 
of  linear  fractional  programming  with  interval  programming  constraints  whose 
matrix  is  of  full  row  rank.  The  analysis  proceeds  by  simple  transformation 
to  canonical  form  exploitation  of  the  Farkas -Minkowski  lemma  and  the  duality 
relationships  which  emerge  from  the  Charnes -Cooper  linear  programming  equi¬ 
valent  for  general  linear  fractional  programming.  The  formulations  as  well 
as  the  proofs  and  the  transformations  provided  by  our  general  linear  frac¬ 
tional  programming  thoery  are  here  employed  to  provide  a  class  of  cases. 

The  argumentation  developing  the  explicit  solution  is  presented,  for  clarity, 
in  an  algor ithmic  format. 


I .  Introduc  tion: 


The  linear  fractional  programming  problem  arises  in  many  contexts  with 
relatively  simple  constraint  sets  --  e.g.,  in  the  reduction  of  integer  pro¬ 
grams  to  knapsack  problems,  in  attrition  games  and  in  Markovian  replacement 
problems  as  well  as  in  Neyman-Pearson  rejection  region  selection  problems. 
Illustrative  examples  are  provided  by  G.  Bradley  [5  ]  or  F.  Glover  and 
R.  E.  Woolsey  [12]— J.  Isbell  and  W.  Marlow  [13],  C.  Derman  [  10]  ,  and 
M.  Klein  [  16]  . 

The  linear  fractional  programming  problem  in  all  generality,  and  with 
all  singular  cases  considered,  was  reduced  in  [8  ]  to  at  most  a  pair  of 
ordinary  linear  programming  problems.  This  immediately  made  available  all 
of  the  algorithms,  interpretations,  etc.,  that  are  associated  with  linear 

2/ 

programming.  This  includes,  we  should  note,  access  to  any  ordered  field,— 
and  any  of  the  algorithms  and  the  computer  codes  for  linear  programming 
problems  which,  by  virtue  of  [8  ]  ,  thereby  also  become  available  for  any 
problem  in  linear  fractional  form.  Thus, with  the  development  in  [8  ]  ,  the 
work  in  linear  fractional  programming  took  a  different  form  from  its  previous 
sole  concern  with  the  development  of  special  types  of  algorithms  for  dealing 
with  this  kind  of  problem. 

In  the  present  paper,  we  apply  our  reduction,  as  given  in  [8  ]  ,  to  a 
general  class  of  linear  fractional  problems  --  viz.,  those  for  which  the  con¬ 
straint  set  is  given  by 

(1.1)  a  <  A  x  <  b 

3  / 

so  that  this  part  of  the  model  is  in  ’’interval  programming"  form.— 

See  also  E.  Balas  and  M.  Padberg  [1  ]. 

See,  e.g.,  the  development  of  the  opposite  sign  theorem  and  related 
developments  in  [7  ]. 

See.  [2  ]  -  [4  ]  and  [18]-[19]. 


3/ 
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Here  we  shall  assume  that  the  matrix  A  is  full  row  rank  and  the  vectors 
a,  b  and  x  meet  the  usual  conditions  for  conformance.  This  means  that  the 
constraint  set  is  a  parallelepiped.  See  the  Final  Appendix  in  [7]. 

Subject  to  conditions  (1.1),  we  wish  to 

T 

(1.2)  maximize  R(x)  =  =  C-,x  jt  Cg  £  constant 

D(x'  dx  +  d 

o 

so  that  we  are  now  concerned  with  a  problem  of  linear  fractional  programming. 
Because  A  is  of  full  row  rank  it  has  a  right  inverse,  A^ ,  and  hence  we  can 
write 

(1.3)  AA#  -  I. 

Now,  setting 

y  =  Ax 


or 


(1.4)  where 


x  =  A#y  +  Pz 


p  =  X  -  A#A 


and 


we  obtain 


z  is  arbitrary 


Tj  T_ 

cAy  +  cPz  +  c 

max  R(y,z)  «  -=rj - £ - ° 

d  Ay  +  d*Pz  +  d 
J  o 


(1.5)  Subject  to 

a  <  y  <  b 

in  place  of  (1.1)  and  (1.2). 

Because  z  is  arbitrary,— ^  unless 

(1.6)  cTP  =  dTP  =  0 


1/ 


Observe  that  we  have  ruled  out  the  case  in  which  R  is  identically 
constant  in  (1.2). 
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we  shall  obtain  max  R  =  °°.  In  order  to  avoid  repetitious  arguments,  however, 

T  T 

we  defer  the  proof  of  this  until  we  have  discussed  the  situation  c  P  =  d  P  *  0. 
See  section  IV,  below. 

Waging  this  consideration  aside,  we  shall  next  proceed  to  solve  this 
problem  explicitly  and  in  all  generality  by  means  of  the  following  3  character¬ 
izations:  The  denominator  D(x)  is  (1)  bisignant,  (2)  unisignant  and  non¬ 
vanishing,  or  (3)  unisignant  and  vanishing  on  the  constraint  set.  In  (1),  i.e., 
the  bisignant  case,  we  shall  show  that  R(x)  =  ®.  Furthermore,  we  shall  show 
how  to  identify  this  case  at  the  outset  so  that  it  may  be  discarded  from 
further  consideration.  This  will  leave  us  with  only  cases  (2)  and  (3)  to 
examine  where  we  shall  proceed  to  transformations  from  which  a  one-pass  numerical 
comparison  of  coefficients  makes  explicit  the  optimal  value  and  solution. 

After  this  has  all  been  done,  we  shall  then  return  to  assumption  (1.6) 
in  a  way  that  utilizes  the  preceding  developments.  Finally  we  shall  supply 
numerical  examples  to  illustrate  some  of  these  situations  and  then  we  shall 
draw  some  conclusions  for  further  research  which  will  return  to  the  remarks 
at  the  opening  of  this  section. 

II ,  bisignant  Denominators: 

Employing  assumption  (1.6),  our  problem  is 


T.#  , 

c  A  y  +  c 


(2.1) 


max  R(y )  =  T  ^ 


d  A  y  +  d 


subject  to 

a  <  y  <  b 

in  place  of  (1.5).  Note,  however,  that  here,  and  in  the  following,  we  shall 
slightly  abuse  notation  by  continuing  to  use  the  symbols  R,  N,  D,  as  in 
(1.1)  and  (1.2)  even  though  we  mean  the  transformed  function,  as  in  (2.1). 

^  A 

Let  D,  D  denote  the  maximum  and  minimum,  respectively,  of  D  over 
the  constraint  set.  We  note: 
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1. omnia  1  : 


(a) 

D 

is 

bisignant  if  and  only  if  D  >  0  and 

D  <  0. 

(b) 

D 

is 

unisignant  if  and  only  if  either  ?) 

<  0  or  D  >  0. 

terms 

of 

y . 

since  we  can  choose  each  component 

of  y  independently 

(1.5) 

-- 

we 

V/  A 

can  express  D  and  D  immediately 

as 

D  =  Z  d  7 ,  b,  +  S  d(.  a .  +  d  =  max  D(y) 

+  j  J  -  J  J  o  VJ,/ 

(2.2) 

D  =  E  d  \  +  Z  d  ^  bj  +  dQ  =  min  D(y) 

where  d'.  s  (dTA^) .  is  the  j—  element  of  dTA^  and  "+"  or  "  indicates 
J  J  J 

that  the  summation  is  on  only  the  positive  or  negative  d  ' j . 

Let  us  first  consider  the  bisignant  situation.  If  we  make  the  trans¬ 
formation  of  variables 


(2.3) 


y.  -  aj  -kjCj,  dj>0 


bj  -  Vj  -  kjCj.  dj<  0 


where  the  k.  >  0  will  be  suitably  chosen,  the  constraints  transform  to 


°  *  kJCj  ^  bj  -  aj 


(2.4)  or 


0  <  Cj  <  6j  a 


b  .  -  a  . 

=  -J _ I 


k. 

J 


Without  loss  of  generality,  the  6^  are  positive,  since  otherwise  Cj  “  0 
and  does  not  enter  into  the  optimization. 


By  choosing  the  k^  suitably,  we  obtain  the  form 


(2.5) 


R(C) 


2  Yj  +  Yo 

J _ _ _ _ _ 

+  j  -  j 


Note  that  E  6^  +  E  6^  >  1 ,  since  otherwise  D  would  not  be  bisignant. 
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One  of  the  following  two  cases  must  now  hold: 

Case  (i)  for  some  £,  0  <  £  <  6,  such  that 
D(£)  =  0  we  have  N(£)  #  0, 
or  else 

Case  (ii)  for  every  £,  0  <  £  <  6,  such  that 
D(£)  -  0  we  have  N(£)  =  0. 

In  case  (i),  since  N(£)  is  continuous,  there  is  a  neighborhood  of  £ 
in  which  N(£)  is  unisignant.  Since 


(2.6) 


Z  £.  +  Z  £.  =  1  <  Z  6.  +  Z  6. 
+  J  -  J  +  J  -  J 


and  0  <  C.  <  6.  in  the  constraint  set,  we  can  choose  e.  >  0,  Z  e.  >  0, 

-  J  -  J  J“  jJ 

so  that  0  <  £  +  e  <  6,  sgn  N(£  +  e)  =  sgn  N(£  -  e) 

and  D(£  +  e)  >  0>  D(£  -  e).  By  approaching  £  along  the  line  segment 

from  one  of  £  +  e.  £  -  e,  we  can  make  R(£)-*®. 

In  case  (ii),  we  must  have  Yj  *  0  for  all  j  such  that  d^  =  0. 

For  D(£)  =  0  involves  specifying  only  the  £^  for  d^  >  0  and  d^  <  0. 

If  Y.  /  0  for  d,  *  0,  then  having  made  D(£)  *  0,  we  can  change  the 
^o  "*o 

value  of  N(C)  by  changing  £.  .  Thus  D(£)  =  0  would  not  imply  N(£)  =  0. 


We  therefore  drop  the  "+", 


1*4.11  II _ II 


notation  in  considering  case  (ii)  and  rewrite 


it  as:  ZyC  +V  =  0  whenever  0  <  £.  <  6.  and  Z  £.  *=  1. 


J  J 


Letting  y^  =  Cjyo»  yQ>  this  becomes, 

i  (?  Yf  y  5  +  Y_  y  )  >  0  whenever 
j  J  J  o  o  - 

Z  y .  -  y  =0 
jf  J  ° 

-  y .  +  6 .  y  >0 
J  J  o  - 


yj 


>  0 


y  >  0 
o  — 


(2.  7) 
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Note  that  tne  implication  extends  to  *  0  since  yo  ■  0  implies  all 


v  =  0 . 
'  ) 


1L 


Wo  now  apply  the  Far kas-Minkowski  lemma-'to  the  pair  of  implications  in 

(2.7)  and  obtain 


.+  .  + 


(2.8) 


f  V.  =  n  -  9.  +  v. 

J  J  J 

1  y  =  +  E  6 .  B*!"  +  v+  ; 


.  e  ,  v  ,  v'  >  o 

J  J  o'  J  J  O  - 


tor  the  first  implication--viz .  ,  (S  y  .  y.  +  y  y  )  >  0.  For  the  second 

J  )  ) 

one  viz . ,  -  (2  y  .  v  .  +  y  y„)  >  0 ,  we  obtain, 

-  j  J  J  o  o 


(2.9) 


- y.  =  pi  -  9 .  -tv. 
J  J  J 


-V  =  -p  +2  6.0.  +  v  ;  0.,  v.,  v  >  0 

o  *  i  J  J  o  j*  j’  o  - 


Adding  the  first  expressions  in  (2.8)  and  (2.9), 


(2.10)  or 


0  =  p+  +  p~  -  (0+  +  Q.)  +  (v+  +  v-) 


6 .  +  0.  **  (u  +  p  )  +  ( V  .  -f  V.) 

J  j  *  ^  J  J 


Adding  the  second  pair 


(2.11)  or 


0  =  -  (p  +  p  )  +  E  6^  ( 9+  +  0^ )  (v+  +  vq) 


p  "bp  -  2  6.  (9.  +  6.)  +  v  +  v 

i  j  J  j  o  o 


Since  each  term  on  the  right  is  non-negative,  we  have 


(2.12) 


p  ■+•  p,  >  0 . 


Next:,  substituting  from  (2.10)  into  (2.11),  we  get 


0  =  -  ( p,  +  p  )  +  E  6.  (p  +  p  )  +  E  6.  (v.  +  v.)  +  v  + 


j  J 


(2.13)  or 


j  '  j  j 


v 


o  o 


0  -  ( Z  6 .  - 1 )  ( p  +  p  )  +  E  6  v  +  E  6 .  v ,  +  v  +  \> 


j  J 


^  u  w  .  V,  I  V  « 

J  J  j  J  j  O  O 


— ^See  Appendix  C  ii  [  7  ] 
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Sincc  the  right-hand  side  is  a  sum  of  non-negative  terms, 
be  zero.  Moreover, 


+  . 

M-  +  M- 


since  E  6.  -  1  >  0 
j  ^ 


(2.14) 


v.  =  v.  -  0  since  6.  >  0 
J  J  J 


v  =  v  =0 
o  o 


By  virtue  of  (2.14),  and  going  back  to  (2.10), 


each  of  these  must 


(2.15) 

Further,  with 

Therefore,  v. 

J 


e+  +  eT  -  o. 

j  j 

9+,  8.  >  0  we  must  have  8"!"  =  0.  ■  0  for  all  j. 

J  J  -  J  J 

*  p,+  ,  for  all  j,  and  yq  =  -p,+  ,  so  that  we  have 


(2.16) 


R(C)  - 


j  J  J 


n+(Z  Q.  -  1) 

j  J 

5£i-x 


+ 


constant. 


In  other  words,  case  (ii)  can  only  occur  in  the  trivial  instance  where  the 
numerator  is  a  constant  multiple  of  the  denominator.  In  this  case,  each 
coefficient  in  the  numerator  is  the  same  multiple  of  the  corresponding 
coefficient  in  the  denominator,  and  this  would  have  to  be  true  in  the  original 
N(x),  D(x)  description  and  hence  obvious  upon  comparing  the  initial  coefficients. 
Since  we  have  ruled  out  this  very  obvious  case--see  (1.2)--we  have  only  max 
R  (x)  =  00  when  D  (x)  is  bisignant  on  the  constraint  set. 

III.  Unisignant  Denominators: 

v 

The  unisignant  cases  now  remain  to  be  considered.  If  D  <  0  we  multiply 
both  N  and  D  by  -1,  (thus  not  altering  the  value  of  R)  and  we  are  then 
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reduced  to  "D"  >  0.  With  this  normalization,  we  make  a  transformation  of 
variables  as  in  (2.3), 


(3.1) 


yj  '  '  8j  ?j  •  dii0 
bj  '  yj  '  gj  -  dj  <  0 


where  the  >  0  will  be  suitably  chosen.  The  constraint  set  will  now  be 


(3.2) 


b  -  a . 

0  <  5.  <  6  -  — 1 - 1 

~  J  “  j  gj 


and,  first  considering  the  case  where 
the  problem  is 


6  >  0, 


the  gj 


can  be  chosen  so  that 


E  Yj  ?j  +  Yo 

(3.3)  max  R  (?)  *  J  E  ?  ~  ’  °  -  -  6j 

j  j 

In  (3.3)  the  summation  is  only  over  "+"  and  "  because,  the  denominator 

being  positive,  optimal  values  for  the  such  that  dj  *  0  can  be 

specified  as  *  0  when  V.  <  0;  5.  -  6.  when  y,  >  0,  and  these  new 

constant  terms  are  assumed  to  be  already  contained  in  y  .  By  the  reduction 

that  we  gave  in  [ff],  however,  the  equivalent  linear  programming  problem  is 

max  E  y.  T1 .  +  y  Tj 
,  J  J  o  'o 


subject  to 

(3.4)  S  T\.  +  *nQ  *  1 


*n.  -  6.T1  <  o 

J  J  O  - 

•Hj  >  o  , 


where  we  can  also  note  that  these  constraints  imply  7]  >  0 

The  dual  to  (3.4)  is 
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min 

subject  to 

(3.5) 


u 


u  +  u).  >  Y- 

j  ~  J 


u  -  2  6.U).  =  Y 
j  J  o 
j 

U)j  >  0 


We  shall  employ  this  dual  in  an  essential  manner  to  obtain  our  desired  one- 
pass  argument  for  obtaining  an  optimum.  At  each  step  in  the  procedure,  we 
shall  have  a  solution  to  a  less  restrictive  problem  than  the  dual  problem 
and  an  associated  primal  feasible  solution. 

Suppose  the  y's  are  renumbered  so  that  Y^  >  Y2  >  •••  Y„ 


n 


Th 


e  case  (i)  yQ  >  has  the  immediately  obvious  primal  solution  Tj^  -  1, 
T|^  *=  0  and  max  R(?)  =  YQ* 

In  the  contrary  case, 

(ii>  V,  >  .  .  .  >  Vp  >  Yc  >  YP+1  >  •  •  •  >  Y„  . 


n 


we  build  up  an  algorithm  based  on  the  dual  problem  in  which  we  choose  u 
t  h 

at  the  q —  step  to  satisfy 


u 


UJ  j  ■  Y  j  .  j  -  1 , .  . .  ,  q 


(3.6) 


uq  -  2  6 .  =  y 

j-i  J  J  0 

Using  the  first  q  equations  to  obtain  in  terms  of  u^  and  substituting 

in  the  last  equation,  we  obtain 


(3.7) 


uM  (1+26 


.)  =  4 


Y. 6.  +  Y 
J  J  c 


and  hence 
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(3.8) 


s  v,6,  +  y 
i  J  J  o 
q  1  J  J 
u  =  - 


q 

1  +  s  6 


j 


Thus,  u"  is  a  convex  combination  of  Yq»  Y^»»«*»  Y^  with  proportionality 

constants  1,  6^,...,  6  .  Note  that  if  uq  <  then  uq,  cuq,  j  «  1 . . 

satisfy  the  first  q  constraints  plus  the  "Y0"  constraint  of  the  dual  problem; 
hence  satisfy  a  less  restrictive  problem  than  the  dual. 

Taking 


(3.9) 


1/(1  +  S  6.) 
1  2 


and 


(3.10) 


T)' 


6./(l  +  Z  6.),  j  -  1 . . 


“Hj  ■  0,  j  >  q  , 

then  T)  j ,  j  *  0,...,  n,  is  a  feasible  solution  to  the  primal  problem  and 

yV  +  YoT)^  =  Uq.  (By  substitution  in  (3.4)  and  comparison  with  (3.8).) 

Hence,  whenever  we  can  get  uq,  uoq  feasible  for  the  dual  problem,  we 
will  have  a  primal  feasible  solution  for  T)q  with  the  same  f unc tional  value 
and  thus  we  will  have  an  optimal  pair  of  dual  solutions .  This,  plus  the 

equivalences  maintained  via  (3.1)  and  our  theory  from  [  ],  thus  justifies 

the  development  that  we  now  detail  as  follows: 


To  start, 
(3.11) 


<Yo  +  Y161)  '  (1  +  6l)  * 


(Note,  u  >  y  since  Yi  >  Y  and  u  is  a  proper  convex  combination 
o  i  o 


of  Y,,  Y0.) 


We  check: 
Is 

If  yes : 


u  >  Y,  ? 


II- 


we  are  done 


and 


If  No.- 


Next  , 


1 

* 

u 

=  u 

1 

* 

U)  . 

=*  uu. 

J 

j 

* 


-  u  , 


> 


j  =  1 
j  >  1 


Tl1  «  1/1+6  x 


T^1-  11 1*  =  6 


1  / 1+6 . 


1  if  9 

ti.  -  *n  *  o,  j  >  i. 


u  <  y2»  then 


2 

u  = 


Vo  +  Vl  +  V2 

1  +  61  +  62 


+  Y,  61 

<■  X  +  6.  >  ( 


1+0, 


1  +  6^  +  &2 


^  +  ^  1  +  6,  +  6  ^  Y2 


1  f  1  +  61 
=  u  ( 


1  +  6^  +  &2 


)  +  ( 


1  + 


61  +  62 


)  v2  <  y2. 


since  u  <  Y2  • 


Is 


u 


>  Y, 


? 


If  Yes:  We  are  done  with  the  substitutions  indicated  by  (3.6)  ff. 


If  No 


2  3 

u  <  and  we  continue  to  u 


This  process  must  stop  by  up  at  the  latest  since 


Y„  +  S  Y  •  6 . 

o  1  j  i 

■p - >  V  >  V  .  > 

1  +  £  6.  o  -  P-11  - 


•  >  X 


n 
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Thus  max  R(^)  =  u  where  s  is  the  least  positive  integer  such  that 


u  >  Ys+l.  and 


u 


s 

V  +  2  Y-  6. 

o  ,  J  j 

s  1  J  J 


s 


1  +  2  6 


j 


This  concludes  the  case  6  >  0. 

The  remaining  case  has  6=0.  By  making  a  transformation  of  variables 


as  in  (3.1),  and  choosing  the  g^  >  0  suitably,  we  obtain 


n 


(3.12) 


f  Yj  +  Yo 

R(5)  -  — - -  .  0  < 


2  ? 


j 


We  may  dispose  of  two  situations  immediately 

(i)  yq  >  0:  then  R(§)  -•  «  as  0. 


(ii) 


YQ  =  0:  here  the  dual  problems  are 


with 


n 


max  2  Y  .  Tj  . 


J  J 


min  u 
wi  th 


=  i 


Tj .  -  6  .  Tj  <  0 

J  J  o  - 


u  +  u).  >  Y- 

J  ~  J 


-  2  6.  0).  >  0 

i  J  J  “ 


V  \  *  0 


u>.  >  0 
J  - 


*  _  _  * 

An  optimal  solution  pair  is  uu^  =  0,  u  =  Y^  =  max  Yj  and  *1^  “  1, 
7)  =0,  j  >  2 .  The  maximum  of  R(§)  is  thus  Y^« 
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The  remaining  instance  is 
(iii) 


Yo<0:  Yl  >.•■>  Yp  >  Yc  >  Yp+l  >.••>  Y„. 


The  dual  linear  programs  can  now  be  written 


max  2  y.  T)  .  4-  y  T) 
i  i  o  o 

j  J 

with  2  "H  .  * 

.  J 


min  u 


with  u  +  u).  >  v. 

J  “  J 


6  ,T)  <0 

J  o  - 


2  6.  u).  =  I  Y 

j  J  J  1  o 


T)j  >  0 


u>.  >  0 

J  “ 


As  before,  we  define 


uq  +  u/? 


V  J 


1 »  •  •  •  >  *1 


1  6.  u/? 

l  J  J 


I  yJ  . 


This  yields 


uq  =  (2  Y.  6  +  Y0)  /  £  6.  . 
1  J  J  °  1  J 


It  may  be  easily  verified  that 


uq  »  u<q'1) 


U  A 

2  6. 


q 

2  6. 

^  V 


( 

q 

2  6. 


(q-1) 


or,  what  is  the  same  thing,  uM  is  a  proper  convex  combination  of  u 

as  in  our  earlier 

r-t-i.  rri 

argument  for  D  >  0. 


,  q  r  _  ,  r  .  r+1  . 

and  u  .  Thus,  if  u  <  Y  ,i»  than  u  <  u  <  y  , 

rtl  r+1 


The  steps  of  our  process  are  as  before:  if  uq  >  Y^+^  we  are  done; 


otherwise,  we  test  u 


q+1 


against  Y^.^*  At  worst  WG  are  done  with 
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n  n  n 

un  =  (E  v  6.  +  V  )  /  Z  6,  . 

1  J  J  °  i-l 

As  before,  uS  =  max  R(5)  where  s  is  the  least  positive  integer  such 
that  us  >  Ys+l  . 

IV.  8  (y.z): 

Returning  to  (1.5)  we  consider  the  remaining  cases  in  which  either 

(a)  drp  =  0,  cTP  *  0 

(4)  or 

(b)  dTP  *  0. 

In  case  (a)  ,  since  z  is  arbitrary  we  can  make  c  P  z  ■  |cPz|-*°°,  hence 

max  R(y,z)  -•  In  case  (b),  we  are  in  the  bisignant  denominator  situation 

T 

since  we  can  make  d  P  z  -*  + 00  .  The  argument  of  the  bisignant  section  of 
this  paper  (with  the  additional  variables,  z)  now  shows  that  max  R  *  ®  since 
we  have  ruled  out,  a  priori,  the  case  in  which  R  s  constant. 

V .  Examples : 

Some  examples  may  help  to  fix  and  sharpen  some  of  the  preceding  develop¬ 
ments.  Thus  consider 

3  -  Xj  +  4 

max  R (x)  =  - — - - 

(5.1)  with 

-1  <  x^  +  X3  <  2 

1  <  x2  <5. 

and  the  variables  x^,  x are  otherwise  unrestricted. 


Here  we  have, 


1/ 


fl  0 


1 


0  0 


(5.2) 


P  =  (I  -  A#  A) 


0  0 
0  0 
_0  0 


1 

0 

0 


To  exhibit  the  development  in  full  detail  we  next  write 


T 

c 


y 


(3,0,  -1) 


1 

0 

0 


cT  P  Z  =  (3,0,  -1) 

0  0-1 

0  0  0 

Z2 

0  0  0 

Zo 

w 

- 

c  =  4 
o 


(0,0, -3) 


(5.3) 

dV'y  *  (0,2,0) 


dT  P  z  =  (0,2,0) 


d  =  0 
o 


1  0 
0  1 
0  0 


f 

* 


— 1 
o 

o 

H  | 

0  0  0 

z2 

0  0  0 

z. 

_  _ 

w 

1/ 


It  may  be  observed  that  the 


choice  is  not  unique. 
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T 

Evidently  in  this  case  d  P  =  0,  as  witness  the  next  to  the  last  expression. 

T  T 

On  the  other  hand,  c  P  ^  0,  as  witness  c  Pz  ■  -  3z^  in  the  second  expression 
for  (5.3).  Hence  condition  (a)  of  the  preceding  section  obtains  and  we  have 
r  -*  oo  even  though 

-i  <  y,  <  2 

(5. A)  1 

1  <  y2  <  5. 


This  occurs  because  z  is  arbitrary  and  can  be  freely  chosen  in 


(5.5) 


TJ  .  T_ 

cAy+cPz+c 


3y, 


max  R  (y ,  z)  = 


3z^  +  A 


d  Ay  +  d  P  z  +  d 


2y. 


which  is  the  specialization  of  (1.5)  to  this  case.  Of  course,  the  result 
r  ->  oo  in  (5.1)  can  be  confirmed  by  direct  inspection,  since  negative  values 
of  x^  may  be  selected  along  with  increasingly  positive  values  of  x^ 
as  required  in  order  to  maintain  the  first  interval  programming  constraint. 
This  last  remark  suggests  that  an  adjunction  such  as 
(5.6)  0  <  x3  <  1 

will  convert  (5.1)  to  a  problem  with  a  finite  maximum.  This  yields  an  A 
for  an  interval  programming  format  with  the  full  row  rank  condition  fulfilled 
as  in 


10  1 

o 

1 

(5.7)  A  = 

0  10 

0  0  1 

0  10 

0  0  1 

_ 

On  the  other  hand,  A  is  also  of  full  column  rank  so  that  we  also  have 

J  A"1  ,, 

A  =  A  and 

p  =  I  -  A^  A  “  0. 

Hence  both  of  the  conditions  specified  in  (1.6)  are  fulfilled  --  viz. 
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T  T 

c  P  =  d  P  =  0 


T  T 

for  any  c  and  d  . 

The  problem  to  be  solved  is  now  written 

T  A# 

c  A  y  +  c 
,  „  o 

max  R  (y)  =  - m - Jl - 

d  A  y  +  d 


3yl  "  4y3  +  4 

2  y. 


(5.8)  with 


-1  <  yx  <  2 
l  <  y2  <  5 
0  <  y3  <  l 


Evidently  the  solution  to  this  problem  is  y^  **  2,  y^  “  0  and  *  1 

so  that 

i  ft 

5 


D  ,  v  10 

max  R  (y)  =  — 


To  obtain  the  corresponding  components  of  x  we  simply  utilize  (1.4)  with 
P  =  0  to  obtain 


4  'N 


(5.9) 


*1 


A  y 


KV 


10-1 
0  10 

0  0  1 


2  ^ 

1 

- 

1 

0  ) 

W 

As  may  be  seen,  these  x  values  satisfy  (5.1)  with  (5.6)  adjoined.  They  are 
evidently  also  maximal  with  R  (x)  =  5  since  x^  can  no  longer  be  negative 
and  X£  and  x^  are  at  their  lower  and  upper  limits,  respectively. 

If  some  cases  the  solutions,  as  above,  may  be  obvious  but,  of  course, 
this  cannot  always  be  expected.  Recourse  to  the  preceding  development,  however, 
will  produce  the  wanted  results  in  any  case,  however,  as  we  illustrate  by 
now  developing  the  above  examples,  along  with  the  related  background  materials, 
in  some  detail  as  follows: 
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itecausc  the  denominator  is  unisignant  we  utilise  section  III.  Observing 
that  d  =  d?  »  2  in  the  denominator  and  hence  is  non-negative* we  have  recourse 
only  to  the  first  part  of  (3.1)  on  page  8  in  order  to  write 


yl  ’  *1  81  h 

(6.1)  y2  *  &2  *  ^2 

y3  -  a3  -  g3  5j 

where,  respectively,  al  -  -1,  a2  -  1  and  a^  -  0^ via  (5.8).  The  development 
from  (3.1)  to  (3.2)  on  page  8  applies  to  this  case  as, 


with 


(6.2) 


0  <  5,  <  6 


.  .b±:...*i  i  -  <-» 


1-1  8,  8, 


b2  "  a2 

o  <  5,  <  62  -  -  ;  *  -  . 
*  1  8*  8 


0  <  §3  <  63 


b,  -  a 


1 


8«! 


The  Insertion  of  (6.1)  into  the  functional  then  produces 


(6.3) 


3  (  81  *  V  *  4  (83  ^3  *  *3>  *  4 

2  (82  ^  *2^ 

3  «i  Si  •  4  «i  ^3  ♦  1 

2  <g2  ?2  +  0 


via  (5.8).  Choosing 


(6.4)  rl  -  1/2,  g2  -  1,  g3  •  1/2 

and  setting 

yx  -  3/2,  v2  -  0,  y3  •  -4/2,  Y0  •  Ul 


(6.5) 
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gives  the  denominator  form  wanted  for  (3.3)  as: 

1.  .i,  +i 

».*  R(?)  . 

«2  +  1 

(6.6)  with  o  <  ?i  <  5 

o  <  52  <  6 
o  <  53  <  2 

The  transformation  5,  “  /T|^  from  our  previously  developed  theorem  [  6] 

then  produces  the  following  example  for  (3.4): 

3  \ 

max  j  0T12  -2  *n3  +  y 

wi  th 

\  +  \  +  ^3  +  \  "  1 

(6.7) 

\  -6T)o  <  0 

■n2  - 0 

*n3  -2*^  <  o 

V  V  ^3  -  0 

where  the  gj  values  of  (6.4)  combine  with  (6.2)  to  give  •  6,  »  4 

and  6^  «  2,  as  required  for  the  application  of  (3.4).  The  corresponding  dual, 

which  our  previous  theory  also  gives  access  to, is 

min  u 


with 

U  +  U), 

>4 

1 

-  2 

(6.8) 

u  + 

U)„ 

>  0 

2 

u  + 

u>. 

>  -2 

3 

u  -6uj^- 

ta2  -  2m3 

uj2, 

~e  ; 

IV 

0. 

This,  of  courset  is  the  application  of  (3.5)  to  the  present  example. 
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since  Yj  ■  2  exceed*  V0  ■  ^  w*  In  the  situation  of  cast  (ll) 
tol lowing  (3.5).  lliu* ,  preserving  our  subscript  identifications'  from  (6.5), 

vo  have 

(6.9)  y{  >Y0>Y2>V3 


in  our  present  situation.  We  therefore  see  that  the  first  application  of  the 
suggested  algorithm  should  suffice.  (See  the  remarks  vhich  conclude  the 
case  6  >  0  in  section  III). 

Applying  (3.11)  now  produces 


(6.10) 


“l  *  (  2  +  2  **)/<!  +  6> 


>  y  m  A 

To  2 


Evidently  this  also  formally  satisfies  the  condition  that  u  equals  or 
exceeds  the  immediate  successor  of  in  (6.9).  Hence,  we  have 


u 


(6.11) 


U) 


u> 


U), 


u 


UL 


IU, 


12 

14 


\  *  U 


3 

2 


12 

14 


JL 

14 


0 

0 


which  satisfy  the  constraints  of  (6.8) ,  as  may  be  verified,  with 
*  . 

min  u  *  u  *  19/14. 


Moving  to  the  primal  problem  via  (3.9), 


1 


1  +  6, 


(6.1  ) 


T) 


1  1+6, 


1  +  6 
6 

1  +  6 


1 

7 

6 

7 


and  all  other  "H  *  0.  See  (3.10). 

Inserting  these  values  for  the  corresponding  T|j  in  (6.7)  we  see  that  all 


constraints  are  satisfied  with 


3  „  'o  3  6.1  1 


19 


a  r\  +  .  a  -  +  —  a  u 

2  1  2  2  *  7  7  *  2  14  * 
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the  same  as  the  value  of  u  ,  thereby  confirming  optimality.  In  fact,  as 
our  theory  [  ]  prescribes,  we  need  merely  apply  the  expressions 

\  =7,  ?!  -  VT|C  *  f  7  *  6 


with  all  other  T|  =  0  and  then  reverse  the  development  from  (6.6)  to  (6.7) 

in  order  to  verify  that  this  value  is  also  optimal  for 

3  ,  1 

2*2  19 

max  R(5)  »  - 6~T“T  “  14 

with 

0  <  t  -  6  <  2  ■<=»  -  6 

"  1  *1 

(6.13) 

0<?3-0<i^-2  -2 

Evidently  we  can  now  directly  effect  substitutions  in  (6.1)  and  obtain 

yl  *  81  ?1  +  al  "  3  "  1  "  2 

(6.14)  y2  *  82  ?2  +  a2  "  0  +  1  “  1 
y3  “  g3  ?3  +  a3  “  0  +  0  8  0 


Then  we  can  proceed  exactly  as  in  (5.9)  to  obtain  the  values  x^  ■  2, 
x2  *  1 ,  x^  ■  0  which  we  previously  observed  to  be  optimal. 
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Conolusion : 

Before  proceeding  any  further  we  should  probably  point  up,  again,  the 
crucial  role  played  by  the  general  theory  (including  the  transformations 
and  proof  procedures)  which  we  introduced  in  [8]  for  making  explicit  contacts 
between  linear  fractional  and  ordinary  linear  programming- -in  all  generality—— 
and  exact  detail.  These  transformations  are  also  utilized  in  the  present 
paper  and  the  theory  is  also  extended  via  the  duality  (and  other)  characteri¬ 
zations  given  in  the  preceding  text.  These  are  joined  together  here  for 
the  proofs  in  algorithmic  format  kind  we  have  just  illustrated  by  example  and 
commentary.  Other  uses  can  undoubtedly  also  be  made  of  this  theory  and  the 
preceding  extensions  via  the  passage  (up  and  back)  between  linear  fractional 
and  ordinary  linear  programming  that  is  now  possible. 

Our  general  theory  has  been  used  by  others,  too,  to  extend  or  simplify 
parts  of  linear  fractional  programming  en  route  to  effecting  the  contacts  with 
ordinary  linear  programming  that  are  thereby  obtained.  The  work  of  Zionts  [22] 
should  perhaps  be  singled  out  as  being  most  immediately  in  line  with  the 
R  =  oo  results  presented  in  this  paper.  Zionts'  development  is  directed  only 
toward  simplifying  matters  by  focusing  on  eliminating  cases  for  linear  frac¬ 
tional  programming  which  are  either  deemed  to  be  unwanted  or  of  little  interest 
for  practical  applications. 

The  developments  cease  as  soon  as  the  contacts  with  ordinary  linear 
programming  are  identified  via  our  theory,  which  he  like  others  utilizes  for 
this  purpose. 

We  have  effected  the  developments  in  this  paper  in  a  way  that  makes 
contact  with  interval  linear  programming.^/  An  opening  for  further  two-way 


i/see  [£),  [3],  [I*]. 
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flows  is  thereby  also  provided.  The  resulting  junctures  should  also  help 
to  guide  subsequent  developments  in  the  more  special  situations  that  now  seem 
to  invite  consideration  in  the  future.  Finally,  the  possibilities  for  dealing 
with  specially  structured  problems  (such  as  those  observed  at  the  start  of 
the  present  paper)  should  also  be  observed  explicitly  in  this  conclusion, 
partly  because  the  theory  we  have  now  developed  and  presented  should  also  be 
a  helpful  guide  to  these  additional  cases  which  are  important  in  their  own 
right.  Thus  we  can  now  conclude  here  by  referring  back  to  our  opening  remarks. 


ADDENDUM 


Although  the  development  in  this  paper  proceeded,  for  clarity,  by 
algorithmic  format,  we  summarize  below  the  explicit  solution  in  tabular 
format  for  direct  theoretical  interpretation  and  utilization. 

Summary  of  Solution 


Transformed 

Numerator 


Bisignant 
Unisignant 
(a)  positive 


(b)  non-negative  y  >  0 

y°  =  0 
To 


s 

u  = 


Y0  +  f  V'2  ,  . 

-  ,  least  s  with  u  ^y  L, 

i  +f  ,2 


y  <  0 
To 


uS  =  Yo  +  T  Y252  ,  least  with  us?y  „ 

1  1  -"TS* .  . .  S+J. 
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